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ABSTRACT 

We argue that bulk spiral flows are ubiquitous in the cool cores (CCs) of clusters and groups of 
galaxies. Such flows are gauged by spiral features in the thermal and chemical properties of the 
intracluster medium, by the multi-phase properties of CCs, and by X-ray edges known as cold fronts. 
We analytically show that observations of piecewise-spiral fronts impose strong constraints on the 
CC, implying the presence of a cold, fast flow, which propagates below a hot, slow inflow, separated 
by a slowly rotating, trailing, quasi spiral, tangential discontinuity surface. This leads to the nearly 



logarithmic spiral pattern, two-phase plasma, p 



density (or T 



temperature) radial profile, 



and ~ 100 kpc size, characteristic of CCs. By advecting heat and mixing the gas, such flows can 
eliminate the cooling problem, provided that a feedback mechanism regulates the flow. In particular, 
we present a quasi-steady-state model for an accretion-quenched, composite flow, in which the fast 
phase is an outflow, regulated by active galactic nucleus bubbles, reproducing the observed low star 
formation rates and explaining some features of bubbles such as their Rb our size. The simplest two- 
component model reproduces several key properties of CCs, so we propose that all such cores harbor 
a spiral flow. Our results can be tested directly in the next few years, for example by ASTRO-H. 
Subject headings: galaxies: clusters: general — hydrodynamics — intergalactic medium — magnetic 
fields — X-rays: galaxies: clusters 



1. INTRODUCTION 

Most galaxy clusters show a central dense, cool core 
(CC) in which the radiative cooling time of the intra- 
cluster medium (ICM) drops well below the age of the 
cluster. Some steady, smoothly distributed heat mech- 
anism is needed in order to balance the observed cool- 
ing, and to sustain the mild (factor of a few, typically) 
temperature drop towards t he center. For rev i ews o f 

Fabian (200G); 



a spiral morphology composed of spatially alternating, 
low entropy and high entropy plasma phases, which are 
approximately at pressure equilibrium. The low entropy 
component is colder and denser (both by a similar factor 
q of up to a few), and higher in metallicity. In some 
cases, the spiral structure may extend beyond the core. 

The Ray leigh- Taylor (RT) stable boundary between 
the low entropy phase from below {i.e. closer to the 
center of the cluster) and high entropy from above forms 
a discontinuity, observed in the form of an X-ray edge 
known as a cold front (CF) when projection effects are 
favorable. Such CFs were found in more than h alf of 
the otherwise relaxed CCs (Markcvit ch et al.ll2003| ). and 
nearly in all of the well-observed, low redshift cores. For 
example, at leas t one CF was found in each of the 10 CCs 
in the sample of Ghizza rdi et al.l ([2010ft . These CFs typi- 
cally show deviations from hydrostatic equilibrium {e.g., 
Mazzotta et al. 2001), a metallicity gradient, and a mild 
pressure jump. In this work we focus exclusively on such 
core-type fronts, and not on the CFs found in merger 
clusters, w here the pressure jump is t ypical ly large (for a 
review, see iMarkevitch fc VikhlinirJ (120071)). nor on the 
putat ive CFs created by shock collisions ([Birnboim et al.1 
12010( 1. where no tangential flow, metallicity, or pressure 
gradients are expected. 

Core CFs are often quasi spiral, piecewise spiral, or 
nearly concentric, suggesting an underlying three di- 
mensional (3D) spiral pattern seen in various projec- 
tions. Multiple, quasi-concentric CFs are often found 
on alternating sides of the center of the cluster, with 
an increasing distance and size, consistent with a spi- 
ral discontinuity manifold observed edge on. Exam- 
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such CC clusters (CCCs), se e Peterson 
iMcNamara fc Nulsenl p007t l: ISokerl (|2010ft 

Energetically, core cooling can plausibly be 
suppressed, for example by thermal conduction 
(Zakamska & Narayan 2003) modified by heat buoyancy 
instabilities (|Quata crt 2008) or by the energy output 
of the active galactic nucleus ( AGN) in the central cD 
galaxy {e.g., IBirzan e"t~a l. 2004). In particular, sufficient 
mechanical en ergy is thought to be deposited in hot 
AGN bubbles dChurazov et all 12002ft. found in at least 
70% of the cool cores (jDunn & Fabianl 12006ft . However, 
it is unclear if the energy can be transferred steadily 
and homogeneously throughout the cooling plasma, in 
order to stem the local thermal instability. Moreover, a 
feedback mechanism is needed in order to regulate the 
heating and adapt to changes in the core. If the AGN 
output is regulated by the accreted, cooling plasma, 
it could furnish such a feedback loop a nd quench the 
globa l thermal instability as well {e.g., IRaffertv et al.l 
12001 . 

In recent years, high resolution X-ray observations 
have revealed the presence of spiral patterns in a 
significant fraction of the CCCs for which high quality 
data ar e available ([Clarke et al.l 120041: iTanaka et al.l 
20061: ISanders et alT~ 

Laeana et al.l 120101: [Randall et al.l 120101: IBlanton et all 
20111: IRoediger" et al j|2011D. These observations suggest 



and A 496 (Markcvitch & Vikhlinin 2007, and references 
therein). The thermal profiles across such CFs indicate 
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that they are strongly shea red tangential disc ontinuities 
(TDs) seen in projection (jKeshet et al.l 120100 . These 
TDs are thought to be stabilized magnetically, with fast, 
nearly sonic, narro w shear layers lying beneath them 
(jKeshet et al.ll20To| ). 

Spiral features and core CFs are typically observed 
only marginally, at low significance, because a high reso- 
lution is needed in order to identify the modest gradients 
and small TD contrasts in projection, and spectral data 
are needed to resolve the temperature/metallicity struc- 
ture. Spirals are rarely unambiguously resolved, but the 
abundance of subtle spiral features in high quality data 
suggests the prevalence of spirals, in particular consider- 
ing projection effects, supporting evidence (such as CFs 
consistent with edge-on spirals) and the complex core 
morphologies (hot bubbles, jets, minor mergers, etc.). 

Moreover, indirect evidence suggests that CCs have a 
universal, unrelaxed structure, which we argue is closely 
related to the spiral features. First, although CCs dif- 
fer in their parameters, reside in different gravitational 
potentials which are strongly affected by the cD galax- 
ies, and are out of equilibrium due to cooling, they 
broadly adhere to universal thermal and chemical pro- 
files. Thus, CCs show radial profiles of density, temper- 
ature and metallicity that cluster aroun d p oc r -1 , T ~ 
r 0A a nd Z ~ r~ 3 , res p ective ly (e.g., iVoigt fc Fabian! 
2001 [Sanderson et al.l 12001 iVikhlinin et al.l 120061: 



Donahue et al.1 2006; Sa nderson et al.ll2009f ). Second, ac- 



cumulating evidence indicates that the plasma in many 
cores, both with and without observed CFs and spiral 
features, cons ists of two (or more) spat i ally separated 
phases (e.g., iMolendi fc Pizzolatoll2001l : iKaastra et al.l 
2004 iWerner et al.ll2010D . These phases are identified as 
associated with the two distinct components that alter- 
nate to give rise to the spiral structure in select nearby 
clusters such as Perseus and Virgo, and can be thus in- 
terpreted in general, as we discuss in [JS] below. 

Core CFs and spiral features were modelled as associ- 
ated with large sca le "sloshing" oscillation s of the ICM, 
driven by mergers t M arkevitch et aLll2001h. possibly in- 
volvin g; only a dark matter subha l o (iTittlev fc Henriksenl 
120051 lAscasibar fc Markevitch] 120061 ). or by weak 
shocks /acoustic waves displacing cold c entral plasma 
(jChurazov et al.l 120031: iFuiita et al.ll2004fl . A subset of 
spiral CFs may be created as low entropy plasma from a 
subcluster core spirals down the gravitati onal potential, 
following an off-axis merger ((Clarke 2004j). 

Sloshing simulations have produced CFs and 
spiral features that n i cely r esemble obser v ations 
(lAscasibar fc Markevitchl 120061: IZuHone et all 120101: 
iRoediger et al.l 120111: iZuHone et al.l 1201 ID . How- 
ever, these features are essentially tr a nsient s 
(|Ascasibar fc Markevitchl 120061 : IZuHone et all I2010I ). 
which do not survive long once co oling and isolation 
by the TD-sheared ma gnetic fields dKeshet et al.l 120101) 
are taken into account fZu Hone et al.l 1201 ID . Moreover, 
these simulations assume a preexisting CC with the 
typically observed size and universal thermal profiles, 
rather than address the formation of such a core. 
Finally, sloshing simulations have not yet addressed the 
role of AGN feedback and its coupling to the flow. 

The preceding discussion suggests that spiral flows are 
ubiquitous in CCs. This is based on the combination of 
the common appearance of spiral structure involving CFs 



and the interpretation of CFs as tangential bulk flows, 
is supported by anecdotal evi dence such as the in tricate 
velocity structure in Perseus (|Sanders et al.ll2004l . figure 
8), and by the agreement between sloshing simulations 
and observed spirals. Moreover, we argue that spiral 
flows may in fact be essential in the formation of a core, in 
shaping its (universal) profile, fixing its size, and avoiding 
catastrophic cooling. Indeed, this would suggest that 
CCs are synonymous with spiral flows. 

Sloshing simulations have demonstrated that a quasi 
stable, spiral mode is easily excited in the core, and de- 
cays only on a long, >Gyr timescale if cooling is not taken 
into account. We analytically study spiral modes in the 
ICM in the presence of cooling and AGN activity. We 
constrain the flow only by the observed, nearly spher- 
ically symmetric pressure profile, and by the presence 
of fas t flows below spiral CFs inferred in iKeshet et al.1 
(2010). The implied features of the core are then shown 
to be strongly constrained, and in good agreement with 
observations and simulations, even for the simplest two- 
component toy model we investigate. 

The paper is organized as follows. In Sj2] we intro- 
duce the framework and the equations governing the flow. 
Power-law scalings of the core and flow parameters are 
derived in [J3J We study the general properties of the 
flow in §31 and explore its two main variants, in which 
the fast flow beneath CFs is either an inflow or an out- 
flow, in [J5l A full solution to one of these variants, which 
combines a hot inflow and a cold, fast outflow, is derived 
in §6l Heating and feedback are discussed, qualitatively, 
in [J7] We examine the observational evidence for a mul- 
tiphase CC plasma in Our results are summarized 
and discussed in Ej9j 

2. DEFINITIONS AND GOVERNING EQUATIONS 
The flow is governed by the continuity equation 

d tP + v-(p\j) = o , (i) 

the momentum equation 

(d t +U-V)U = a = -/;- 1 VP + g , (2) 
and the energy equation 



PT 



ds 



dt T - 1 



(dt + XJ- V)T + pTV-U 

= V-(kVT) - A 



(3) 



Here, we denote the mass density p, the pressure P, the 
temperature T, the normalized temperature T = kgT/p, 
the Boltzmann constant ks 1 the mean particle mass 
p ~ 0.6m p , the proton mass m p , the velocity vector U, 
the inertial acceleration vector a, the entropy s, and the 
thermal conduction coefficient k. 
For simplicity, the cooling function A is approximated 

as 

< A oc p 2 T 1/2 , (4) 

independent of metallicity; the effect of a different tem- 
perature dependence is discussed at the end of fj3] We 
approximate the equation of state as that of an ideal gas, 



pk B T 
P = = pT , 



(5) 
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with adiabatic index T = 5/3. For simplicity, we assume 
that the gravitational potential $ is spherically symmet- 
ric, such that the gravitational acceleration g = — V$ 
is radial. In addition, we neglect thermal conduction by 
taking k = 0. This is probably justified for conduction 
perpendicular to the flow, as the shear amplified mag- 
netic fields quench the perpendicular transport. Parallel 
condu ction is suppressed by sca ttering off plasma waves 
{e.g., Pistinncr & Eichler 1998); the role of parallel con- 
duction is discussed in SJ7] 

The observation of a spiral pattern in a galaxy cluster 
suggests a symmetry axis, denoted by z, perpendicular 
to the plane of the spiral. We assume that the spiral flow 
evolves sufficiently slowly to approach a steady state, in 
a (possibly rotating) cluster frame of reference. There- 
fore, if the spiral pattern is rotating, it must do so uni- 
formly, in order to avoid a winding problem. Henceforth, 
we assume that the pattern is rotating with a constant 
but arbitrary angular frequency u about the z axis, and 
work in the corresponding rotating frame unless other- 
wise stated. This requires the addition of fictitious forces 
to the momentum equation, which then becomes 



in a simple form: continuity, 



(U- V)U = 



-p _1 VP - 
-2ft x U 



- ft x (SI x x) 



(6) 



where ft = ujz is the rotation vector. 

For simplicity, we assume that an equatorial plane can 
be found, perpendicular to the z symmetry axis, in which 
all streamlines are approximately confined to the plane. 
We use cylindrical coordinates x = {r, <j>, z}, with r being 
the radial distance and <f> the azimuth, defined such that 
the origin coincides with the center of the cluster, and 
the equatorial plane lies at z = 0. We label the velocity 
components by U = {v, w,U 2 }, and sometimes use the 
two-dimensional (2D) velocity parallel to the equatorial 
plane, u = {v,w}. The TD pattern in such a plane per- 
pendicular to z is quantified by the angle a(r) it makes 

with respect to <fi. These definitions are illustrated in 
Figure [TJ 




Fig. 1. — Parameter illustration in a plane parallel to the equa- 
tor. The TD (solid curve) forms an angle a with cf>, i.e. with the 
circle r = constant (dot-dashed). The plasma must flow parallel 
to the TD. The planar velocity u (thick green arrow; parallel to 
the TD at the point analyzed) is decomposed into radial (v) and 
angular (w) components (dashed green arrows). The pattern may 
rotate about the z-axis, in the +z sense if lo > (circular arrow). 



The flow along a TD must be parallel to it. Thus, 
focusing on the flow along (on either side of) a TD in the 
equatorial plane, the steady state equations can be cast 



A, 



din p 



d\nr v 
momentum conservation along 

d I w 

dr \ 2 J r 7 
momentum conservation along r, 



V-U = - 1 + A W + A, 



2 ''■ W 2 O0 



2ww ; 



v 

dr \~2 

and energy conservation, 
vpT 



= a r + u r + 2ww ; 



(7) 



(8) 



(9) 



1 + A w + A T + A A = -A . 



(10) 

Here, r = T 1 '^^ 1 ' is the reduced temperature, 7 = 
sin (a) serves as the pattern parameter, full derivatives 
are taken along the flow, 

d_ + J_d_ 

dr jr d(f> ' 



d 
dr 



(11) 



and we defined power-law indices 

X A , ^ (12, 
dmr 

for the evolution of each quantity A along the flow. These 
equations apply along streamlines that are parallel to the 
TD, both above and below it, everywhere that Eq. (jlip 
holds. Deviations are expected far from the TD, but for a 
spiral flow these are significant only in a localized mixing 
layer interleaved between the TDs. 

In order to incorporate the 3D structure of the flow in 
a simple manner, we assume that U and the TD surface 
are perpendicular when projected onto the z — r plane. 
Parameterizing the effects of the flow component per- 
pendicular to the plane of interest is thus simplified, as 
it depends only on the TD radius of curvature in the r — z 
plane, Re = -[1 + r'(z) 2 ] 3 / 2 /r"(z), through the relation 



1 v) 77T 



(13) 



This assumption is illustrated in Figure [5J 

Equivalently, this assumption means that the flow is 
confined to the plane spanned by the local TD normal, 
n, and n x z. This generalizes the equatorial analy- 
sis which follows to any point in the core, both inside 
and outside the equatorial plane, as long as streamlines 
within the confined plane remain parallel to the TD. Note 
that a flow of this type is well-defined globally only if 
r/Rg < 1, such that the flow planes do not intersect. 
Moreover, a pure equatorial spiral pattern is in general 
deformed in such confined planes that are not parallel to 
the equator, as 7(7") may become complicated and even 
non-monotonic. 

Next, we solve the tangential momentum equation (JHJ, 
by assuming that the inertial acceleration is approxi- 
mately radial, \a$\ -C \a r \. This is justified by the ap- 
proximately spherical distributions of pressure and of to- 
tal (including dark matter) gravitating mass in CCs, al- 
though we shall later see that first order corrections in 



4 



U. Keshet 




Fig. 2. — Illustration of the flow perpendicular to the equatorial 
plane, and of CF projection along the line of sight i (black arrow). 
Three examples of CFs, intersecting the equatorial plane (dotted 
horizontal line) at the same distance r from the center (black disk), 
are shown in a plane perpendicular to the equator, i.e. "edge on". 
The CFs have radii of curvature Rg = r (solid red), Rg > r (dashed 
blue), and Rg < r (dot-dashed green). The locations (open circles) 
where the three CFs are parallel to i and thus most easily observed 
in projection, lie at separations (double arrow lines) respectively 
equal to, longer than, and shorter than 2r. The flow is assumed, for 
simplicity, to be confined to planes (long-dashed red lines shown 
for the case Rg = r) spanned by the local TD normal, n, and n X z, 
where z is the symmetry axis (vertical arrow). 



|a0/a r | must be retained. For a 
two solutions: 

1. The trivial solution, 

v = w = 



0, Eq. 



admits 



(14) 



describes corotation with the pattern. In the iner- 
tial (non-rotating) frame, this becomes w"' = tor. 



2. A non-trivial solution, 



w = 



■ tor ■ 
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C w 
r 



LOT 



(15) 



where C w is a constant. Here, in the inertial frame, 
W W — C w /r. For example, for C w = 0, the inertial 
frame motion is purely radial: the gas propagates 
inward/outward radially along the rotating spiral, 
in resemblance of an Archimedes screw. 

3. SCALING OF A TWO-COMPONENT FLOW 

An important clue to the nature of core spirals is 
the presence of fast, nearly sonic flows, found beneath 
well-observed CFs. Such flows, and the shear they 
induce across and below the TDs, are needed (see 
iKeshet et alj |2010) in order to explain the thermal pro- 
files observed across CFs, the spatial correlation be- 
tween CFs and radio miniha los ([Mazzotta & Giac intuccl 
120081: IKeshet fc Loebl 120101) . and the remarkable sta- 
bility of the TDs, witnessed by the ubiquity, smooth- 
ness, and sub-mean free pa th thinness of the CFs 
(|Markevitch & Vikhlinml 120071 ). 

3.1. Nearly adiabatic fast component 

Although the cooling time is shorter than the age of the 
cluster, it is much longer than the radial sound-crossing 
time, even deep in the core. Therefore, the fast flow is 



approximately adiabatic. Equations (JSJ), (J7J, and (|10[) 
then yield for this fast (indices /) component: 



X,™ = ^ = ^ - -0.36 , 
At (/) = Ll± X = -A ~ -0.24 



w + Ro 
where we took 



r + A 



A = A, 



-0.64 



-0.6 



(16) 



(17) 



(18) 



(19) 



on the right hand side (RHS) of these equations as a 
typical value of the logarithmic pressure gradient. To il- 
lustrate the weak dependence upon A, Table[T]reproduces 
these and the following relations for a slightly different 
pressure profile, A = —2/3. 

Next consider the momentum equation along </>. The 
trivial solution Eq. (fT4| corotates with the TD pattern, 
so the fast component must follow the other, non-trivial 
solution, Eq. (TT51) . Moreover, if the wr term in Eq. (|15[) 
were to dominate the fast component anywhere, then 
Eq. (ITSl would yield Rg < 0, corresponding to outward- 
curved CFs, unlike the nearly concentric CFs typically 
observed. We conclude that the flow beneath the TD 
must be faster than the pattern throughout the flow, such 
that Eq. {T5J yields 

A w ,/ -1 , (20) 

as expected from angular momentum conservation. 

Eq. (|18[) now indicates that the perpendicular radius 
of curvature is larger than r, 



r 

Re 



A 

T 



4a 

5 



0.36 



(21) 



so the TDs are intermediate between spherical and cylin- 
drical. Namely, for an edge-on (i.e. within the equato- 
rial plane) observer, the CFs would appear intermediate 
between concentric semicircles, and lines parallel to the 
z-axis. Such an intermediate curvature is needed in or- 
der to explain the observation of nearly concentric CFs, 
arguably seen in various projections. It is also needed for 
the flow parameterization in Eq. (fT5|) to be well-defined 
outside the equatorial plane. Note that a similar per- 
pendicular TD curvature , of order r/Rg ~ 1/2, is found 
in sloshing simulati ons (jAs casibar & Marke vitchl 120061 : 
IRoediger etl^l201lT ). 

3.2. Radial momentum conservation poorly constrains 
the power-law model 

The momentum equation in the r direction can be in- 
tegrated to yield a Bernoulli equation, 



f~l T/ 



7 



(22) 



where C$(0) is an arbitrary function, u = |u|, and 
we neglected the w-dependent terms by assuming that 
M ^ l w l/ r - However, this equation is not directly use- 
ful for deriving the power-law scaling relations. As the 
flow is assumed subsonic, the second term on the left 
hand side (LHS) is subdominant (by at least a factor 
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of three). As $ rises monotonically with r throughout 
the core's gravitational potential well, Eq. ([22]) is locally 
consistent with a radially declining temperature profile 
of the fast flow component, as derived above. It becomes 
increasingly difficult to reproduce such a declining tem- 
perature profile away from the core, where the gravita- 
tional well becomes shallow: this may confine such fast 
spiral flows to the vicinity of the core. 

The slow part of the flow, which strongly dominates 
the core, is nearly in hydrostatic equilibrium. For a hy- 
drostatic plasma with a temperature profile T oc r A r^o.4^ 
Eq. ([22]) is thus satisfied for the fast component ifl 



AT , 



(23) 



up to an additive constant. As C$ should be independent 
of r, this calibration holds only locally. However, the 
additive constant and the weak radial dependence of T on 
the RHS of Eq. (|23|) indicate that the radial momentum 
equation can be satisfied globally, by slightly relaxing 
the approximations leading to Eq. (f22|) . including the 
velocity term in Eq. (|23[) . or allowing deviations from a 
pure power-law behavior. 

In principle, Eq. (|2"2"| can be solved for the spiral pat- 
tern 7(r), as 



Wy \ I 



1^ 



(24) 



However, the 7 term is a small correction here (7 2 ~ 0.04; 
see i}6]). Therefore, in addition to the reservations men- 
tioned above, the resulting pattern would be sensitive 
to the precise form of $, would depend on the cj correc- 
tions, and in general would not produce a pure power-law 
profile. Instead, we shall compute ~f{r) using the slow 
component of the flow. 

3.3. Nearly corotating slow inflow 

Consider the slow flow just above the CF. If this com- 
ponent were to follow the non-trivial solution Eq. (|15[) . 
then there would be fast flows above the CF, at ei- 
ther small or large radii. Such flows are not observed; 
rather, the region above the C F is typically con s istent 
with hydrostatic equilibrium (Markcvi tch et al.1 120011 : 
iKeshet et alJ l2010). We deduce that the velocity of the 
slow component is approximately given by the trivial so- 
lution Eq. (TT4"]), If the flow were adiabatic, this would 
correspond to stationary, hydrostatic equilibrium plasma 
in the rotating frame, with an arbitrary thermal profile 
and consistent with any spiral pattern. 

However, a slow flow, in particular one with v ~ 0, 
must be modified by cooling. Indeed, a stationary flow 
with u = does not satisfy the energy equation (|TU)) . 
Therefore, we must consider small deviations from a triv- 
ial flow and from a purely radial acceleration. Linearizing 
Eq. (jHJ to first order in {a0,v,w} yields 



2ujv ~ 



1 



(25) 



—d^P , 
rp 

1 One may alternatively model the gravitating mass dist ribution, 
for example using an NFW profile (Navarro et al. 1996) for the 
dark matter, but the cD galaxy must be included due to its strong 
influence on the core. 



so it is natural to assume the scaling 

X v {s) ~ A — 1 - A/ s ) , (26) 

where indices s pertain to the slow component. Note that 
the linearized Eq. ([9]) then implies a small deviation from 

hydrostatic equilibrium, ai^ ~ —u) 2 r + (p s r r y)~ 1 d c j > P. 

Finally, combining the continuity equation (|7|) , the tan- 
gential momentum equation (|2"6"|) . and the energy equa- 
tion for the cooling, slow component, 



A v « + 1a t W _ l * A P W , 



(27) 



along with the perpendicular radius of curvature deter- 
mined by the fast component in Eq. (|21j) and the equa- 
tion of state (|5|), fixes the scaling of the slow flow: 



A (s) 



-4 + 3A 



1.16 



A T « * * 0.56 



and 



A 



(») 



5 

-1 + 2A 



-0.44 



(all three independent of T), and 

1 + 3A A _ 1 
5 + f ~~ ~5 



A (s) 



(independent of A for T — 5/3). 
These relations also fix the pattern of the spiral, 



A 7 = A V ( S >-A W (^^— -A = ^A 
1 5 



-0.24 



(28) 
(29) 

(30) 
(31) 

(32) 



The pattern in turn may be used to determine the scal- 
ing of the radial velocity of the fast flow, which was left 
undetermined above: 



a v (/)^^a-i = Ha-i 

1 5 



-1.24 



(33) 



These power-law indices are all fairly weak functions 
of A, indicating that the model is robust. They are sum- 
marized in Table [TJ evaluated numerically for a slightly 
different pressure profile, A = —2/3. 



TABLE 1 

Scaling of the two-component model, for r = |; A : 



Property A 


Power law index A^ = din \A /dln(r) 


pressure P 
spiral slope 7 


A = Ap = — | 
§A=-^~-0.27 




Fast phase Slow phase 


density p 
temperature T 
rad. velocity v 
tan. velocity w 


j?A = I = 0.4 -*+ 3A = 1= 1.2 
|A = -^~-0.27 ^# = -,1^0.53 

5 15 5 15 

|A 1- r? ~ 1.27 -!+ 2A - 7 ~ 0.47 

5 15 5 15 

-1 -i = -0.2 



a Note that A here is taken slightly different numerically than in 
the text, in order to show its rather weak effect, in particular on the 
slow component. The flow equations require r/Rg = — A/T = 2/5, 
intermediate between cylindrical and spherical TD geometry. 

The above scaling relations imply that the bracketed 
term in the energy equation (fT0|) is positive for the slow 
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component, 

1 + A w (s) + A« + X z ~ [4T + (5 - 3T)A]/[5(r - 1)] ~ 2 . 

Therefore, the slow component must be an inflow, 
with v s < 0. Approximating the cooling rate 
as A ~ 2.1 x lQ-^nl cm-^TlK] 1 / 2 erg s^ 1 cm- 3 
(jRvbicki fc Lightmanl lT9861 we find 

v s ~ -10 n om T^ 1/2 r 10 km s" 1 , (34) 

where no. 03 is the electron number density n = p/fj, 
in 0.03 cm" 3 units, T 4 = [k B T/A keV), and r w = 
(r/10 kpc). This velocity is of order 1% of the sound 
velocity c s , consistent with the slow flows above CFs in- 
ferred from observations. 

Notice that the spiral pattern and fast flow are inde- 
pendent of, and the slow flow depends weakly on, the 
precise temperature dependence of the cooling function 
Eq. For example, in the cold, high metallicity re- 
gions in the centers of small clusters and groups of galax- 
ies, t he temperature dependen ce weakens to A oc T° 
(e.g., iPeterson fc Fabianl |2006[ ) . Here, the scaling laws 

in Eqs. ([25 ]) -([3"I ]) slightly change, to X p {s} ~ -1.07, 
X T {s) ^ 0.47, A v (s) ~ -0.53, and X w {s) ~ -0.29. 

4. GENERAL PROPERTIES OF THE FLOW 

The scaling relations derived in JJ] indicate that the 
velocity of the fast, cold component increases at smaller 
radii, while the slow, hot flow becomes denser and some- 
what cooler at small radii. The spiral pattern is found to 
be nearly logarithmic (for a logarithmic spiral 7 = const., 
or equivalently, A 7 = 0), with a slight tendency towards 
an Archimedes spiral (for which A 7 = — 1) . We may now 
combine these scaling laws with observations in order to 
test the model and furth er explore the natu re of the flow. 

Observations suggest dKeshet et al.ll20iol) that the fast 
flow inward of a CF is confined to the vicinity of the 
discontinuity, inducing shear up to distances of order 
Ar < 0.2r below the CF at radius r. Therefore, in 
the simple approximation where the core is composed 
of only the two, fast and slow, components, we expect 
the thermal properties of the core to be dominated by 
the slow flow. Hence, the core must approximately fol- 
low p ~ r" 1 and T ~ r 1 / 2 . Similarly, the behavior of the 
TD shear is dominated by the fast flow, so approximately 

Alt = U s — Uf r~j 7'- 1 . 

These properties are broadly consiste nt with the ther- 
mal st ructure of observ e d core s (e.g.. iVoigt fe Fabianl 
2001 [Sanderson et all I200H IVikhlinin et all 1200a 
Donahue et all 120061: iSanderson et all 120091). and with 
the ( presently weak) constraints on shear (jKeshet et"all 
2010). We also find that the spiral pattern inferred 
above is consistent with the observed spiral pattern in 
Perseus, which is amongst the clearest spirals observed 
to date, and is presumably seen nearly face-on (approx- 
imately along the z-axis); see Figures [3] and 0] Even 
better agreement is found with the late-time, quasi- 
stable pattern discovere d in t he merger simulation of 
lAscasibar fc M arkevitch (2006), where the projection is 
controlled; see Figures [5] and [5] In both cases, the spiral 
is approximately logarithmic, with slightly negative A 7 . 

The two components of the flow slide across each other 
along the spiral, RT stable, TD surface that separates be- 



tween the fast, cold flow from below and the slow, hot 
component from above. This strongly sheared boundary 
is thought to be stabilized magne tically, isolating the two 
components (jKeshet et al.l 12010). In contrast, the tran- 
sition layer between slow flow from below and fast flow 
from above is not discontinuous, as such a discontinuity 
would become RT unstable. This layer is likely to involve 
substantial heat conduction, convection, and gas mixing, 
and may become turbulent. In this region our analysis 
breaks down, as streamlines may no longer be parallel to 
the TD. 

In order for the TD to remain RT stable, the fast flow 
below the discontinuity must remain at least as dense 
and as cool as the flow above the TD. The temperature 
(or equivalently, as we assume an isobaric transition, the 
density) contrast q across the TD increases with r as 

K^jpL^ocr 4 / 5 , (35) 
T f p s 

so the discontinuity would become RT unstable at very 
small radii, limiting the extent of the flow or the validity 
of the scaling derived above. Similarly, the velocity ratio 
between the fast and slow components diminishes with 
increasing radius, 

1 < — ~ =F— oc r~ 4 / 5 , (36) 
u s u s 

where the minus (plus) sign corresponds to a fast out- 
flow (inflow). Thus, the shear responsible for the stabi- 
lizing magnetic fields gradually diminishes at large radii, 
again limiting the extent of the flow or the validity of the 
power-law scaling. Note that the ±4/5 power-law indices 
in Eqs. (|35[) and (|36| are independent of T and A. 

As TD contrasts of order q ~ 2 and nearly sonic shear 
Ait ~ c s are inferred from CF observations, extrapolat- 
ing these scalings to small radii suggests that the two 
components have similar temperatures and densities at 
some small radius ro, near the base of the TD, where q 
approaches unity and the fast component is nearly sonic, 

q(r ) ~ 1 ; it/(r ) a c s . (37) 

More generally, a measured contrast q(r) or shear 
Ait(r) thus places an upper limit on the radial extent 
of the discontinuity, for example ro > g(r)~ 5 / 4 r. As an 
illustration, a TD following the power-law scalings of $3] 
extends no more than a factor r/ro ~ 2.4 (~ 3.9) below a 
CF observed at radius r, if the inferred TD contrast there 
is q(r) = 2 (= 3). In practice, the flow may extend some- 
what beyond such estimates, because (i)the flow can 
persist beyond the TD, if the discontinuity evolves into a 
more gradual transition; (ii) inaccurate deprojection may 
cause q and Au to appear smaller than they really are, 
for example if a very dense or fast flow is confined to the 
very near vicinity of the TD; and (Hi) deviations from 
our power-law model are expected at small radii, due to 
the nearly sonic velocity of the fast component, the pos- 
sible breakdown of the assumption u s = |u s | <C u>r, and 
the active environment near the AGN. 

The scaling relations derived in iJ3] are likely to break 
down at some large radius r max , where the fast flow be- 
comes sufficiently slow that its cooling can no longer be 
neglected. We may use Eqs. ([3lJ). ([36]). and ([3Tj). to 
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Fig. 3. — Spiral pattern observed in Perseus. Shown is an XMM image (Churazov ct al. 2003, data courtesy: E. Churazov), before (left; 
30' diameter) and after (right; 25' diameter) smoothing with an adaptive, photon count per-bin conserving, Gaussian filter, cropping edge 
effects, and applying a gradient filter. Red contours enclose dark regions that correspond to enhanced gradients, tracing CFs. The light 
regions have small gradients, and show mosaic and other artifacts. The dashed green curve is the best-fit power-law spiral derived in Figure 
□ 

{e.g., Richtmyer-Meshkov) instabilities, turbulence, and 
the weaker seed magnetic fields anticipated at large radii, 
this suggests that the two components mix at some radius 
r > r max - Both effects may contribute to the absence of 
observed CFs with large contrasts. 

The momentum equations for the slow component re- 
quire that to =fi 0, i.e. the spiral pattern must rotate, 
albeit much slower than the fast component. The sense 
of rotation is related to the tangential force acting on the 
slow component, by Eq. (with v < 0): if points 
towards the outspiral (inspiral) direction, then the spiral 
is trailing (leading). For a spherically symmetric gravita- 
tional potential, sign(a^) = — sign(d,pP), so the sense of 
rotation is fixed by the minute tangential pressure gra- 
dients observed. 

In a quasi spherical system, the deviation of the pres- 
sure map from its azimuthal average may reveal such 
small gradients. In Perseus, such a procedure shows a 
slightly elevated thermal pressure significantly outside 
the CFs, at a finite distance above them: compare the 
pressure map in figure 10 of IChurazov"eT"aLl (pOOl to 
the temperature or entropy maps. This corresponds to 
a pressure gradient within the slow component point- 
ing towards the inspiral direction. It suggests that the 
spiral pattern is (%) indeed rotating; and (ii) trailing, be- 
cause the force ~ exerted on the slow component is 
directed along the outspiral direction. The above con- 
clusions are not sensitive to the presence of a magnetic 
pressure, because the elevated thermal pressure is found 
significantly far above the CF: one does not expect the 
magnetic pressure to increase to significant levels as one 
approaches the CF from above. 

The projected appearance of a spiral TD depends on 
the orientation of the spiral axis with respect to the line 
of sight, and on the 3D structure of the TD manifold. As 
CFs are often nearly concentric, and some spiral CFs are 
observed, Re/r cannot be much larger or much smaller 
than unity. Some projection effects are illustrated in Fig- 
ure [5] (a full analysis of projection is deferred to future 
work). The figure shows how, for Rg > r (in dashed 
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Fig. 4. — Discontinuity pattern in Perseus. The strong gradient 
regions (enclosed by red contours in Figure |3)l are shown in the r—0 
phase space. Although 8 is 27r-periodic (pale regions), one may 
trace it along an extended interval (dark regions) by following the 
pattern in Figure [3] and varying the threshold. The best fit (red 
curve) power-law spiral, 6>[rad] ~ —107. 8 + 78. 3(r/arcsec) 054 , cor- 
responds to A 7 = —1 — r9"(r)/9'(r) ~ —0.05: nearly a logarithmic 
spiral. 

crudely estimate that 

W ^ (S) V2 * 100 (—) V2 k P c • (38) 

The spiral flow may persist beyond r max with a some- 
what different scaling, or may retain the same scaling if 
the cooling of the cold component is compensated, for 
example by heat conduction from the hot component. 
Note that by r max , the contrast across the TD becomes 
large, of order q ~ (r max 

Alternatively, the TD may be disrupted by Kelvin- 
Helmholtz (KH) and other instabilities near or somewhat 
above r max . At small radii, shear magnetic amplification 
is thought to generate s trong magnetic fiel ds that sta- 
bilize the discontinuity (|Keshet et al.ll2010l ). At larger 
radii, the shear gradually weakens. This implies weaker 
KH instabilities, on one hand, but slower magnetic am- 
plification, on the other. Taking into account additional 
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blue), as the line of sight gradually deviates from face- 
on viewing (in which the line of sight is parallel to the 
z axis of symmetry), the spiral pattern is stretched 'par- 
allel to the projected z axis. An opposite, perpendicu- 
lar stretching occurs for Rg < r (shown in dot-dashed 
green). Due to the poor statistics, it is difficult to mea- 
sure this curvature observationally and test Eq. (f2"Tj) and 
the corresponding parallel stretching. 

We have determined the radius of curvature of the flow 
perpendicular to the equatorial plane, but did not ad- 
dress the evolution of streamlines far from the equator, 
nor have we specified the extent of the spiral flow above 
the plane. The flow must extend a substantial fraction 
of the way towards the poles, in order to account for the 
observation of spiral patterns at various projections. In 
principle, the flow can extend all the way to the pole, 
where the model streamlines have an approximate heli- 
cal pattern around the z-axis, with an amplitude that 
diminishes close to the pole. 

5. TWO TYPES OF COMPOSITE FLOWS 

Due to the adiabatic nature of the fast component, 
its behavior is symmetric under time-reversal. Hence, in 
contrast to the slow component which must be an inflovfi 
the analysis above does not distinguish whether the fast 
component is an inflow or an outflow. We thus consider 
below both options for the composite spiral flow: (i) a 
combination of slow and fast inflows; or combina- 
tion of a fast outflow and a slow inflow. 

The former type of flow is essentially a transient that 
accretes mass onto the center, while the latter could in 
principle be long-lived, and may form a closed circulation 
loop with quenched accretion. Note that although spi- 
ral inflow and outflow components carry opposite signs 
of mass and heat radial fluxes, they carry the same sign 
of angular momentum flux. Hence, both an inflow and 
an outflow similarly transfer angular momentum in the 
trailing (leading) spiral sense inward (outward), across 
the core, limiting the lifetime of the flow unless an exter- 
nal torque is involved. 

5.1. fast and slow inspiral 

The first option, in which the fast component is an 
inflow, is in essence a transient phenomenon (albeit long- 
lasting, of timescale > Gyr in the absence of cooling) , as 
the cold, metal-rich inflow must originate deep in the 
core. In order to produce such a flow, the central plasma 
must first be uplifted by some violent process, such as 
the gravitational kick and a ngular momentum depos ition 
imparted by a merger event (iMarkevi tch et al.ll200lD. the 
impact of weak shocks o r acoustic waves (jChurazov et al.1 
12001 iFujita eTall200l . or the uplift of hot (radio) AGN 
bubbles that buoyantly rise from the center. 

In particular, sloshing motions induced by strong 
perturbations of the core following a merger event 
were shown to produce CFs and spiral patterns in 
broad agreement with observations, provided that the 
merger does not strongly disrupt the original core. 
This was demonstrated in simulat ions of major merg- 
ers with a dark matter subhalo (|Tittlev &: Henriksenl 

2 This is due to cooling. A slow outflow is possible, for example, 
in an adiabatic simulation. 
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Fig. 5. — Spiral flow formed as a result of a mass ratio 5 merger 
wi th a dark matter subclu s ter, in the adiabatic SPH simulation 
of Ascasibar & Markevitch (2006, data courtesy: Y. Ascasibar). 
The slice shown is 560 kpc on the side. Greyscale shows the gas 
temperature, ranging from 3.3 keV (dark) to 9.8 keV (bright). Ar- 
rows trace the gas flow. Green contours enclose regions in which 
(V X U) z is smaller than a negative threshold, i.e. strong local 
clockwise circulation. The dark matter peak is shown as a yellow 
triangle. The spiral pattern rotates counterclockwise, i.e. is trail- 
ing. The fast component here is an inflow. The slow component 
is in part an outflow, although this is obscured by the uniform 
rotation of the pattern. 




Fig. 6. — Discon t inuity pattern in the simulation of 
Ascasibar & Markevitch (2006) shown in Figure [5] derived in the 
same method as Figure [4] based on regions exceeding a gradient 
threshold (blue points). The best power-law spiral fit (red curve), 
6 = -26.2 + 8.5(r/kpc) ' 27 , corresponds to A 7 ~ -0.27: nearly a 
logarithmic (slightly tending towards an Archimedes) spiral. 



l2005tlAscasibar fc M arkevitch 2006|) and of minor merg- 
ers (jRoediger et al. || 2011|). and in simulatio ns that incor- 
porated cooling (jZuHone et al.ll2010L 1201 lh . 

However, without additional effects, such sloshing can- 
not protect CCs from catastrophic cooling. Although 
shown to somewhat delay cooling by mixing the gas 
(|ZuHone e t al. 2010), sloshing alone cannot sustain the 
core for more than ~ 1 Gyr, once magnetic shear ampli- 
fication along the TPs ( Keshet et al.ll2010T ) is taken into 
account (jZuHone et al.l 20111 ). Nevertheless, it remains 
to be seen if AGN feedback could stem cooling in the 
presence of such a composite spiral inflow. 

Strong magnetic fields, which are thought to form 
not only at the TD, but throughout the fast flow layer 
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llKeshet et al.l I2010T) as r ecently demonstrated numeri- 
callv (jZuHone et al ll20lTI ). can divert AGN bubbles and 
entrain them within the spiral pattern. Such entrain- 
ment is indeed observed in som e cores, as illustrate d 
for Perseus i n Figu re [7J See also iForman et al.l (|2007| ): 
iFabian et al.1 {2011). Projection effects confuse the in- 
terpretation of such images, in particular in regards to 
the interaction of the spiral flow with entrained bubbles, 
linear cold filaments, and other localized features. For 
example, some bubbles axes or long filaments may in 
fact be oriented perpendicular to the equatorial plane, 
along the z symmetry axis. 

In addition to enhancing the dissipation of bubble en- 
ergy, such interaction between rising bubbles and a spiral 
inflow suggests a natural feedback mechanism: accretion- 
driven AGN activity may disrupt or heat the inflow. An- 
other possibility is that mechanical feedback may stop 
the inflow and push it outwards, resulting in an inflow- 
outflow composite spiral of the type studied in ij5.2l be- 
low. Feedback is further discussed in SJ7J 

5.2. fast outspiral and slow inspiral 

The second possible type of composite flow combines a 
fast, cold outflow and a slow, hot inflow, spiraling across 
each other. Such a flow follows the same scaling relations 
as the double-inflow model of £15.11 although shear and 
gas mixing are probably somewhat stronger. 

It is interesting to note that adiabatic sloshing simu- 
lations sometimes show, at least locally, composite flows 
of an opposite character, combining a fast inflow and 
a slow outflow, as seen for example in som e regio ns in 
the snapshot of the lAscasibar fc Markevitchl (|2006t ) sim- 
ulation depicted in Figure [5] The flow discussed here 
is similar, in the adiabatic limit, to the time-reversal of 
such an adiabatic simulation. 

A novel property of a composite inspiral-outspiral flow 
of the type discussed here is that, unlike the composite 
inflow model of £15.11 it could constitute a quasi steady- 
state solution, as long as the center can absorb angular 
momentum in the trailing spiral sense. This can occur, 
for example, if an external force torques the center in the 
opposite, leading spiral direction. Feedback-regulated 
flows that alternate between fast and slow mass accre- 
tion, or between spirals of an opposite orientation, are 
also possible. 

Specifically, consider the mass flux across a given ra- 
dius in the equatorial plane. The outflow layer is much 
faster than the inflow, but is also much narrower, so the 
relative widths of the two components (defined, say, as 
the ratio v(r) = A(f>f/A(j) s < 1 of their extents at a 
given r) can be chosen to eliminate the radial mass flux. 
Due to the curvature of the streamlines both inside and 
outside the equatorial plane, the planar mass flux pvr of 
each component is r-dependent. Nevertheless, for both 
components it follows the same scaling, 



dln(pvr) r 



(39) 



Thus, a fixed ratio v between inflow and outflow vol- 
ume fractions would guarantee a choked mass accretion 
throughout the flow, if v is properly tuned. Note that, as 
the streamlines are approximately parallel, such a spiral 
flow carries no angular momentum in the rotating refer- 



ence frame, although it does transfer angular momentum 
radially, across the flow region. 

Such a non-accreting flow essentially constitutes a 
closed circulation loop, which could expla in the low star 
formation rate observed in cores (e.g., Raffertv et al.l 
[200l IQ'Dea et al.l[2008l) . Here, cold gas from the base 
of the flow at ro is rapidly removed from the center out- 
wards, decelerates, and in part slowly accretes back to 
the center, after being non-adiabatically heated by close 
proximity to hotter gas, through gas mixing and heat 
conduction. The flow thus solves the cooling problem, 
by advecting heat from the huge heat reservoir outside 
the core and removing dense gas from the center before it 
can catastrophically cool, while quenching or maintain- 
ing a low level of mass accretion, and carrying angular 
momentum away from the center. 

Such a flow requires some mechanism to rapidly propel 
gas from the central parts of the core outwards. Can- 
didate mechanisms include gravitational forces, such as 
the stirring effect of the rotating vector separating be- 
tween gas and dark matter peaks, and the AGN energy 
output, coupled for example through the forces exerted 
by shocks, jets, and buoyant hot bubbles. Such AGN 
bubbles, in particular, suggest a natural feedback mech- 
anism, provided that their production is moderated by 
the accretion rate, as discussed in iJTJ 

High angular momentum in the inner core is another 
requirement, needed both in order to sustain the an- 
gular momentum flux, and to loosen the gravitational 
binding of the core. Such angular momentum can build 
up through (possibly alternating) epochs of strong spi- 
ral flows, for example in the form of an inspiral-inspiral 
composite flow, by merger events, and by the sinking 
of dark matter substructure towards the center. The 
initial formation of the spiral pattern can be facilitated 
by a merger event, by amplified g-modes, or by spiral 
instabilities that tend to separate even a non-rotating, 
uniformly accreting, cooling gas into spiral motions of 
opposite orientation, for example in the standing ac- 
cretion shock ins t ability (SASI) simulations discussed in 
IBlondin fc Shawl $MY7\ ). 

6. TWO-COMPONENT INFLOW-OUTFLOW TOY MODEL 

Above we derived some of the properties and scaling 
relations of a spiral flow of the type expected in a cluster 
core, constrained by the presence of a spiral TD surface 
with strong shear due to a fast flow below the disconti- 
nuity. We showed that the flows immediately above and 
immediately below the TD follow different scaling rela- 
tions, but have not quantified the variation in the flow 
properties with distance from the TD. Here we wish to 
provide the full solution to one type of such a spiral flow, 
by introducing a simple two-component toy model. 

In order to obtain a fully constrained system and re- 
produce the low mass accretion rate inferred from obser- 
vations, we choose the flow variant of ij5.21 combining a 
fast outflow and a slow inflow, and require that the ra- 
dial mass flux vanishes. In order to account for the mild, 
q ~ 2 contrast of typical CFs, and the nearly sonic shear 
inferred from their profiles, we assume that pf = p s , 
Tf = T s , and Uf = |u/| = c s at some small radius ro, as 
discussed in SJU see Eq. (p?7|) . These constraints fix all the 
model parameters except the normalization of 7, which 
can be approximately measured as demonstrated below. 
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Fig. 7. — Deep (900 ks ACIS-S) archival Chandra image of Perseus (presented, e.g., by Fabian ct al. 2006; Sanders & Fabian 2007), 8' 
(~ 240 kpc) on the side, shown with (right) and without (left) the features discussed in the text. Over-plotted on the right are the best-fit 
spiral CF pattern (dashed blue) derived in Figure [4] along with apparent (green closed curves) and possible (yellow curves) X-ray cavities, 
suggesting AGN bubbles interacting, and possibly entrained in the spiral flow. The X-ray c avities were identifie d by a visual inspection 
(M. Markevitch, private communications, 2009, 2011); some of them were also pointed out in Fabian ct al. 1 2 01lT) . 



The requirement of choked mass accretion fixes the ra- 
tio between the widths of the two layers, which we pa- 
rameterize using the angular extents Acf> of each phase 
at a given r, 



function of the normalized radius a = r/ro by 



p(a) oc A<f>f(j 



A m 



(42) 



oc 



A0/ 
A0 S 



c s 7 



0.05 



/ «0.03?To 



V T470.2 

"(40) 

Here we defined 70.2 = 7( r o)/0.2 ~ 1, based on pre- 
liminary estimates of 7 at the base of the spiral pat- 
terns: (i) in Perseu s, where 7(10 kpc) ~ 0.20 (see F igure 
[3]); and (ii)'m the lAscasibar fc M arkevitch (200fj|) sim- 
ulation, where 7(10 kpc) ~ 0.23 (see Figure [3]). Recall 
that although v is evaluated at rg, it remains constant 
along the spiral; cf. Eq. (|39|) . 

A logarithmic spiral follows r(0) oc e 7< ^. As the spi- 
ral patterns observed are nearly logarithmic, the radial 
separation between consecutive TD spiral arms in the 
equatorial plane is approximately r n+ i/r n ~ e 2,7 * 1 ~ 
3.5e 27r ' 7-0 ' 2 '. In the toy model, the fast flow is confined 
to a thin layer below the TD, with an approximate radial 
thickness of 



4/5 



0-05i/o.o5 ( 



1 



4/5 



-1.16 



where ^0.05 = (^/0.05), and we used a pressure profile 
A ~ —0.6 in the last line. This profile is best fit for 
the above typical parameters by p oc cr -11 in the range 
1 < g < 10, as shown in Figure [8] (left panel). 

The emissivity-weighted temperature profile is simi- 
larly given by 



T(cr) oc 



A<f)f<J 



|A T <«+2A <« 



A(j> s (T 



Aififa 



A(j) s a 



iA r < = >+2A>) 



1 - V (1 - a 2 / 5 ) 4+2A 

° C 1 + ^(ct6/5-1) CT 

^ 1-0.05^0.05 (1-g 2 / 5 ) Q 

~ 1 + 0.05^0.05 (<J 6/5 ~ 1) a 



(43) 



5d 



Ar f ~ r (l 



-27T7iA 



0.06r 



(41) 



where we used the parameters of Eq. ([4"U|) . In reality, 
the fast layer may be more extended, and the transition 
between the fast and slow phases must be gradual. 

Next, we compute the azimuthal averages of the flow 
properties within the equatorial plane. Note that al- 
though these quantities approximately give the projected 
radial profiles when the streamlines lie on spheres [i.e. 
when R.g ~ r), a viewing angle-dependent correction is 
needed to account for projection effects when Rg ^ r, in 
particular when the spiral is not observed face-on. Nev- 
ertheless, as long as Rg/r is not too far removed from 
unity, such corrections are in general small. 

The azimuthally averaged density profile is given as a 



The logarithmic slope of T(a) declines with a. The best 
fit power law for the above typical parameters is T oc a 0A 
in the range 1 < a < 5, as shown in Figure [8] (right 
panel). The figure also shows, for comparison, the uni- 
yersal mean deproject ed temperature profile found by 
ISanderson et al.1 (|2006[ ) . scaled with ro = 0.01r 50 o, where 
rsoo is the radius enclosing 500 times the critical density 
of the Universe. 

These density and temperature profiles are in good 
agreement with the characteristic structure of cores. In 
particular, the temperature peak typically observed near 
''max is reproduced, although the power-law solution may 
break down above ~ r max due to cooling and the steep- 
ening of the pressure profile. The azimuthally averaged 
profiles depend on the volume fraction v of the fast flow: 
smaller (larger) values of v lead to a steeper (flatter) 
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Fig. 8. — Equatorial radial profiles of density (left panel) and temperature (right panel) in the accretion-choked spiral flow model 
described in §E\ Shown are the azimuthally averaged profiles (solid curves), their best-fit power-laws p oc r~ 11 and T oc r 0A (dot-dashed 
curves) evaluated near rg, and the profiles of the dominant, slow (short dashed) and subdominant, fast (long dashed) co mpon ents. Th e 
model parameters used here are no. 03»"lo/T4 = 1 and 7(ro) = 0.2, or equivalently u = 0.05 and r max /ro = 10; see Eqs. J38D and (1406 . 
The power- law solution is probably modified beyond ~ r ma x (dotted horizontal line), unless cooling of the fast component is balanced, for 
example by heat conduction. We assume that T = 5/3 and A = —0.6; i n reality the pressure gradient A steepens above ~ r ma x- The mean 
deprojected CCC temperature profile found by Sanderson ct al. (2006) is shown for rg = O.Olrsoo (dotted curve in the right panel, with a 
band representing the la confidence level). 



density decline and temperature rise, approaching (devi- 
ating from) the slow phase profiles shown as black short- 
dashed curves in Figure This may explain the scatter 
seen around the universal thermal profiles among differ- 
ent clusters. 

Various properties of the model can be directly read off 
Eqs. (1551) - (|5T|) . and (fltfl) . in particular within the equa- 
torial plane. For example, the fast outflow contains a 
small, v ~ 5% fraction of the mass at the base of the 
spiral r , but an increasing, [1 + (y~ x — l)(j"/Vo) _4 ' 5 ] _1 
fraction at larger radii. The kinetic energy of this fast 
outflow is a small, - (Tv/S)(u f /c s f ~ 0.03(r/r ) _8/5 , 
radially decreasing fraction of the thermal energy. 

Eq. (|40j) . which guarantees a vanishing mass accretion, 
similarly implies that the spiral flow carries no angular 
momentum in the rotating frame, because the slow and 
fast streamlines are parallel. The angular momentum of 
the flow region is therefore L ~ ojI, where the moment 
of inertia is I ~ (7r/8)r 2 M(< r) in the p oc r _1 approxi- 
mation. 

As mentioned above, although spiral flows may carry 
negligible angular momentum, for example in the choked- 
accretion case with lu — > 0, the flow necessarily transfers 
angular momentum radially, because (i)the angular mo- 
mentum fluxes of the two components always add up, 
rather than cancel; and (ii) here the angular momen- 
tum flux carried by the fast component strongly dom- 
inates that carried by the slow component, by a factor 

Uf/u a ~ 20(^/0.05)- 1 CAo)~ 4/5 - 

Thus, angular momentum in the sense of a leading spi- 
ral rotation is transferred outwards, within the equatorial 
plane and perpendicular to it, towards the poles. Equiva- 
lently, angular momentum in the sense of a trailing spiral 
rotation is transported inward. This torques the central 
region in the trailing direction, for example slowing it 
down if it initially rotates in a leading sense. This an- 
gular momentum inflow sets a characteristic timescale 
during which the equatorial rotation velocity at the base 
of the spiral changes by the sound speed, 



\VfWfVj V T 470.2/, 



y 0.05 



Gyr . (44) 



A significant fraction of the angular momentum inflow 
may be transferred to the uniform rotation of the spiral 
flow region, thus extending t^, as demonstrated in |J7]for 
the case of torqued bubbles. However, as the uniform 
rotation must remain slow, the spiral flow must eventu- 
ally stop or reverse, unless angular momentum is lost to 
the cluster's periphery or an additional source of angular 
momentum torques the center in the opposite, leading 



7. THERMAL AND MECHANICAL FEEDBACK 

An energy source that offsets cooling in the core, and 
a feedback loop that regulates it, are thought to be es- 
sential in quenching the global thermal instability. AGN 
output and heat conduction can in principle stem cool- 
ing, provi ded that the feed back efficiency is sufficiently 
high (e.g., iGuo et a l. 2008), and that the energy is prop- 
erly distributed in space and time. Here we discuss the 
nature of feedback in general terms; a quantitative anal- 
ysis is beyond the scope of the present study. 

Heat conduction and convection can be significantly 
elevated in the spiral flow region, due to the influence of 
the bulk flow on the transport processes and the close 
proximity between hot and cold phases. Magnetic shear 
amplification leads to strong magnetic fields parallel to 
the streamlines, along the TD and within the fast shear 
layer beneath it, thus suppressing tran sverse transport 
(jKeshet et alj|2010t IZuHone et al.ll2011l) . However, par- 
allel transport may be enhanced due to the stronger, 
aligned fields, with less plasma wave inhibition. More 
importantly, the RT unstable layer is not expected to 
undergo coherent shear; it could support substantial ra- 
dial heat conduction and convection. 

In order to close the feedback loop, the AGN output 
must be regulated by mass accretion through the slow 
inspiral, as well as by the fast component if it is an 
inflow (in the composite inflow case discussed in §5.1|) . 
Although the fast flow is denser, it is not sensitive to 
heating (or cooling) due to its fast, near-sonic velocity, 
in particular at small radii. In contrast, the cooling reg- 
ulated, slow component is susceptible to heating, but the 
resulting behavior is sensitive to the details of the heat 
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deposition. Importantly, both fast and slow components 
can be significantly modified mechanically, for example 
by the forces exerted by shocks, jets, and in particular 
buoyant AGN bubbles. 

Consider heating first. Our model remains valid, with 
catastrophic cooling being quenched, if heating of the 
slow inflow is both efficient and able to approximately re- 
tain the X v ^ ~ — 1/2 scaling derived above. This may be 
achieved, for example, if the specific heating rate roughly 
scales, on average, as e(r) ~ r~ 2 , like the A oc p 2 T x l 2 
cooling profile. These requirements are less strict than 
the e(r) oc p 2 scaling sometimes invoked, as they rely on 
heat conduction and gas mixing for distributing the heat 
along streamlines and on small scales. 

Next, consider a mechanical forcing of the flow. Slow- 
ing the inflow would only have a cooling effect, as it would 
increase the ratio between the inspiral time and the cool- 
ing time. In the limit where the forcing is strong enough 
to cause the infall to change sign and become an out- 
flow, but still remain slow and susceptible to cooling, a 
steady flow would imply unrealistic thermal profiles, for 
example A p > A T (c/. Eqs.lTland lTOl) . 

We conclude that mechanical forces can effectively 
stem cooling in a (quasi) steady-state flow of the type 
studied here, only by pushing gas outwards at high ve- 
locity. The fast, outgoing flow may then drag some of 
the inspiraling gas along, causing the remaining inflow 
to expand (thus slightly cooling adiabatically), and the 
cooling rate to decline. The ultimate outcome, even for 
an initially pure inflow, is a slow inspiral-fast outspiral, 
composite flow, of the type discussed in i j5.2l and ^Hl with 
an outflow-to-inflow volume ratio v regulated by the me- 
chanical feedback. 

In this scenario, the inflow-outflow spiral composite 
flow becomes a direct manifestation of cooling and feed- 
back, in the presence of angular momentum. The slow 
inflow, which dominates the core, is in essence a cooling 
flow, whereas the fast outflow removes dense gas from 
the center of the core before it has a chance to catas- 
trophically cool. Both components and the flux ratios 
between them are regulated by feedback, self-consistently 
enforced by the spiral-constrained flow. The steady state 
here corresponds to a closed circulation with quenched 
accretion, advecting heat inward and angular momen- 
tum (in the leading spiral sense) outward. The small, 
few percent ratio between the kinetic energy of the out- 
flow and the thermal energy (see S}6]) indicates that this 
mechanism is energetically plausible. 

The rise of AGN bubbles is a natural candidate for me- 
diating the necessary mechanical feedback, as the bub- 
bles preferentially deposit their energy at or just below 
the strongly magnetized, cool layer. Preliminary evi- 
dence for this is seen in the suggested entrainment of 
bubbles within the spiral pattern; see Figure [71 

The presence of the spiral, with its ordered bulk flow, 
shear, and magnetic layer, constrains the motion of bub- 
bles through the flow, and facilitates the transfer of ki- 
netic bubble energy to the inflowing gas. In particu- 
lar, the spiral directs the bubbles, by channelling them 
along the magnetic layer and furnishing them with angu- 
lar momentum in the leading spiral sense. Thus, the gas 
encountered by the bubbles acquires more momentum 
along the outspiral direction, and less random motions, 



than it would obtain in the absence of a spiral. This is 
clear in the extreme case where bubbles are sufficiently 
large to span a substantial fraction of the cross section 
area of the inflow, such that gas is either forced outwards 
or compressed below the bubbles. This would agree with 
the large bubble size Rb and i ts tight, linear scal ing with 
distance, Rb ~ 0.6r, found bv lDiehl et al.l (|2008l see fig- 
ure 1); note that the inflow cross section is indeed similar 
in size and in radial dependence, ~ e _27r7 r/(l + j/) ~ 0.7r 
(for a logarithmic spiral; see SJSJ) . Conversely, confine- 
ment by the spiral magnetic structure can explain the 
observed size and size-distance relation of bubbles. 

This confinement, so far ignored in feedback studies, 
is likely to play a major role in the flow. The force that 
the spiral structure exerts on the bubbles, and the com- 
pression of gas below the bubbles, could explain the en- 
hanced pressure observed above the TD, and the deduced 
trailing rotation of the spiral (see It may also ex- 
plain the high angular momentum of the outflow near 
the base of the spiral: bubbles are pushed in the leading 
direction by the magnetic layer already near or below ro, 
and throughout the spiral, while the more massive spi- 
ral structure thus acquires equal angular momentum in 
the trailing rotation sense. Thus, once the magnetic spi- 
ral structure has formed, e.g., by a SASI-like instability, 
the flow can in principle be sustained (well beyond it 
in Eq. |44|) by mechanical feedback without an additional 
source of angular momentum, and in particular without 
invoking mergers, if the trailing rotation imparted to the 
spiral structure is transferred to large radii. Such a sce- 
nario, in which there is neither a net angular momentum 
nor a net angular momentum flux, is beyond the scope of 
the present analysis (e.g., it probably breaks the uniform 
rotation assumption), and will be examined numerically 
in future work. 

8. OBSERVATIONAL EVIDENCE FOR TWO COMPONENTS 

Observations indicate that CCs typically harbor more 
than one plasma phase at any given distance from the 
center of the cluster. This is seen quite generally, in- 
cluding in clusters where no spiral patt ern was so far 
detected (see, e.g., iKaastra et al.l 12004) . The charac- 
teristic tem perature spread at a given radius is a factor 
of - 2 - 4 (|Buote et all 120031; IKaastra et all 120041 ). see 
Figure [9] (abscissa) . An identical spread exists in den- 
sity, as azimuthal pressure gradients are negligible. A 
single phase is typically observed only in the outer parts 
of the core, although this may be due to poor statistics 
(|Kaastra et al.ll2004f) . 

The azimuthally averaged temperature, and the tem- 
perature of the hottest phase, were shown to decline 
towards the center of the cluster, approximately as a 
power law, T oc t Xt with \ T ~ 0.4 (jKaastra et al.l 120041 : 
Piff aretti fc Kaastrall2006l). 

iMolendi fc Pizzolatol (|2001l ) claim that X-ray observa- 
tions of CCs are consistent at any given radius with a 
single plasma phase, with azimuthal temperature varia- 
tions. There have been attempts to attribute the multi- 
ple temperature measurements in some clusters to pro- 
jection effects through an in herently spherical t empera- 
ture profile (e.g., in A478, Ide Plaa et all 120041) . How- 
ever, this was shown not to be the case in select cases, for 
example in 2 A 0335+096 (jTanaka et al.ll2006l ). Indeed, 
in some well-resolved clusters, coherent azimuthal tern- 
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perature variations of a factor 2-4 are clearly observed, 
sometimes taking the form of a spiral, as discussed in SJTJ 
The emission measure was found to scale roughly as 
Y = j p 2 dV ~ T aY with an average ay ~ 4, in various 
clusters and different radii (Kaastr a et al.l [2004) . The 
hotter phase or phases thus strongly dominate the distri- 
bution. For example, for two phases at pressure balance 
with an average temperature ratio (q) = T max /T m - m ~ 2, 
this would imply a small volume fraction of the cold 
phase, 



a Y +2 



0.02 



-6- (ay -4) 



(45) 



where the result strongly depends upon the poorly con- 
strained (q) and ay- 

We point out in Figure El that the results of multiphase 
modelling (jKaastra et al.ll2004T ) suggest a correlation be- 
tween the range T m i„/T max of temperatures observed at 
a given radius, and the spatial variation Ti n /T out of the 
temperature of the hottest phase across the core. (How- 
ever, the statistical errors are substantial.) Studies of 
individual clus ters also suggest th at T m ; n ~ Ti„ (e.g., in 
2 A 0335+096, iTanaka et al1l2006D . Such a correlation is 
consistent with a model in which the ICM at any radius 
in the core combines hot plasma from outside the core, 
and cool plasma of a nearly constant temperature. 




Fig. 9. — Ratio between the inner and outer annul us temperature 
of the hottest phase in CCs ( data from tables 6- 7 inlKaastra et all 
1(2001 ): Perseus data from ISanders & Fabianl (20071 )). plotted 
(green diamonds) against the average ratio T m ; n /Tmax between 
the coldest and hottest ph ases in the same ann ulus (average of 
best fit column in table 8 of Kaastra et al. (2004)). Higher quality 
data are highlighted (blue circles with linear fit in red, for clusters 
with data for > 3 annuli). Statistical errors (not shown) are large; 
we use only data with standard deviation smaller than half the 
mean. 

While the hottest, dominant phase usually becomes 
colder towards the center of the cluster, the behavior 
of the cold, subdominant component/s is less clear. In 
some clusters, the data are consistent with a cold phase 
of constant temperature. Thi s is illustrated for V irgo 
in Figure [TU1 using data from Kaastra et al.l (|2004 we 
chose Virgo because its bins are the most numerou s, and 
amongst the least noisy) and I Werner et all (|2010f ). 

The above lines of evidence, combined at face value, 
suggest azimuthal temperature variations by a factor of 
2-4, throughout the core, which is on average strongly 



dominated by a hot, T a 



J3.4 



phase, with a small (of 



order a few percent by volume) contribution from a cold, 
T m in ~ const, component. These are indeed the qualita- 
tive features derived in the two-component spiral model. 



^" 2 - 
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Fig. f 0. — Projected radial temperature profiles of the hot (filled 
symbols) and cold (empty symbols) phases in Virgo. Shown are the 
azi muthally averag e d tem peratures (circles; data from tables 6 and 
8 of Kaastra et al. (2004)), and the temper atures measur e d alon g 
the southwestern X-ray arm (squares; from Werner et al. (2010)). 
Here, 1' corresponds to ~ 4.7 kpc. 

Note, however, that the significance of the evidence re- 
viewed in this section is in general quite low. 

9. SUMMARY AND DISCUSSION 

We argue that spiral flows are ubiquitously present 
in cool cores, and show that the properties of the CCs 
and of the spiral flows are strongly constrained by the 
very appearance of piecewise spiral cold fronts in high- 
resolution X-ray maps. We interpret such core CFs as 
projected tangential discontinuities that lie above layers 
of fast flowin g, cold gas, as infer red from the observed 
CF profiles in lKeshet etail (|2010D . 

In order to obtain a simple, analytic model, we make 
several simplifying assumptions/approximations: (i) a 
symmetry axis z, with a preferred (equatorial, spiral) 
perpendicular plane; (ii) a steady state flow in some ro- 
tating (about z) frame; (Hi) an ideal gas of adiabatic in- 
dex r = 5/3; (iv) a spherically symmetric gravitational 
potential, i.e. in particular, negligible self gravity of the 
spiral pattern; (v) an approximately thermal pressure, 
nearly (but not exactly) spherically symmetric, scaling 
as P cx r A --°- 6 (our results depend weakly on A); (vi) a 
two-component flow approximation; (vii) negligible heat 
conduction and convection; (viii) Row outside the equa- 
torial plane is confined to the plane spanned by the TD 
normal h and n x z, so the equatorial analysis trivially 
generalizes to the entire volume; and (ix) d<j,P ~ r x . 

Under these assumptions and approximations, we find 
that the TDs and streamlines must follow a bulging cylin- 
drical spiral manifold: parallel to the equatorial plane, 
the manifold is approximately a logarithmic (tending 
slightly towards an Archimedes, A 7 < 0) spiral, while 
perpendicular to the equator it is intermediate between 
a concentric semicircle and a line parallel to the z axis; 
see Figures Q] and [2] Pictorially, one may image such a 
TD surface as a (rather useless) open barrel, in which 
the metal hoops were replaced by spiral bands. The TD 
parameters (curvature Re perpendicular to the equator, 
and spiral pattern A 7 ) are summarized in Table [T] The 
TD pattern must slowly and uniformly rotate around 
the z-axis. Subtle pressure gradients (observed at a fi- 
nite distance above CFs in Perseus) suggest that the spi- 
ral is trailing. The spiral pattern remains apparent in 
obliquely projected CFs, although it is increasingly dis- 
torted and broken to pieces as the viewing axis deviated 
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from z (see Figure [2]). 

We examine the spiral pattern observed in Perseus (see 
Figures [3] and |U) and the pattern found in the slosh- 
ing simulation of lAscasibar fc Markevitchl (|2006l gener- 
ated by a major merger with a dark subcluster; see Fig- 
ures [5] and O • Indeed, in both cases a nearly logarith- 
mic spiral pattern is found. The pattern in the simu- 
lations, in particular, is an approximately uniformly ro- 
tating, trailing, logarithmic (slightly tending towards an 
Archimedes) spiral, in agreement with the model. 

In order to account for CF observations, the spiral flow 
must combine a dense, cold, and metal-rich component, 
flowing fast and nearly adiabatically in a narrow layer 
beneath the TD, and a slowly inspiraling, cooling, hotter 
component, which fills most of the volume. These two 
flow components correspond to the two distinct solutions 
of the angular momentum equation (jSJ. Extrapolating 
the properties of observed CFs suggest that at some small 
radius ro, the two components have comparable densi- 
ties and temperatures, and the fast component is nearly 
sonic, such that q ~ (r/r ) 4 / 5 and uj/u s ~ 20(r/ro) _4//5 , 
independent of Y and A; see Eqs. (|3S])-(|37]) and gO]). A 
transition layer that separates between the two compo- 
nents at their RT unstable side (where the fast, dense 
flow lies above the slow flow) mixes the gas and may 
heat the cold component. 

The scaling properties derived for each component of 
the flow are summarized in Table [TJ The density and 
temperature of the fast, adiabatic component change 
slowly with radius, and it is increasingly fast near the 
center. In contrast, the velocity of the slow component 
changes more gradually, and it becomes much denser and 
somewhat colder at small radii. The azimuthally aver- 
aged density and temperature profiles are dominated by 
the slow component, and are found to be in good agree- 
ment with observation^: p ~ r" - 1 and T ~ r Ar ~ 05 . The 
model qualitatively agrees with the multiphase proper- 
ties observed in CCs (see Sj8|), and reproduces the typical 
~ 100 kpc core size; see Eq. (|38|) . 

The equations governing the flow do not break a degen- 
eracy in the direction in which the fast, adiabatic compo- 
nent is flowing, allowing for two types of composite spi- 
ral flows. In the first type (see §5.1|) . both fast and slow 
components are inflows, such as expected in the slosh- 
ing oscillations induced after a merger event disrupts the 
core (in the presence of cooling). This is essentially a 
transient mode, with substantial accretion and little re- 
silience to cooling, unless AGN feedback reverses part of 
the slow infall and accelerates it outwards to high veloc- 
ity. In such a scenario, the core evolves into a second 
type of flow: an inspiral-outspiral composite mode. 

This second type of composite flow (see £ !5.2j) involves 
a cold, fast outflow and a hot, slow inflow. While this 
requires some mechanism for launching or sustaining the 
outflow, it naturally solves the global cooling problem 
with a low level of mass accretion, provided that the 
outflow is regulated by some accretion-driven feedback 
mechanism. AGN bubbles, in particular, provide both a 
launching mechanism and a plausible feedback loop; see 
the discussions in Sj6j SJ7J and below. 

3 One of the two independent thermal parameters is essentially 
fixed by the gravitational potential; the second is a prediction of 
the model. Here we choose to externally fix the pressure profile. 



An additional assumption of choked mass accretion 
may account for the weak cooling flows and the low star 
formation rates inferred from observations. This selects 
the second, inflow-outflow type of composite spiral flow, 
strongly constrains the properties of such a flow, and 
produces thermal profiles in good agreement with obser- 
vations (see Figure [5]). Such a flow is essentially a closed 
circulation, carrying no angular momentum in the rotat- 
ing frame, advecting heat inward from outside the core, 
and carrying cold, dense gas outwards at high velocity 
before it can catastrophically cool. 

Both types of spiral flows transfer angular momentum 
radially, torquing the central core in the sense of a trailing 
spiral rotation. This can persist over timescales longer 
than tL ~ 1 Gyr (see Eq. ((44])) only if the spiral pattern 
reverses, additional angular momentum is deposited in 
the center, or the angular momentum involved in the 
trailing rotation of the entire structure is transferred 
{e.g., magnetically; see outwards. 

Both types of flows may quench the local thermal in- 
stability, as they efficiently distribute energy spatially 
and temporally across the core by bringing cold gas from 
small radii and hot gas from outside the core into close 
contact, strongly mix the gas in the RT unstable layer, 
and possibly enhance transport along streamlines. 

Moreover, in the presence of AGN feedback, the spiral 
may stabilize the core also against the global thermal in- 
stability, as discussed in fjjj One way this could happen 
is through heating, but this requires an approximately 
r~ 2 heating profile. A more natural and robust feed- 
back involves the mechanical energy output of the AGN, 
forcing gas outwards at high velocity. This can lead to 
stable flows of the second, inflow-outflow type, with re- 
alistic core profiles and choked mass accretion. Here, the 
spiral plays a central role in coupling the AGN output 
to the gas, for example by channelling bubbles along the 
spiral. This efficiently quenches cooling by regulating the 
inflow, removing dense gas from the base of the spiral, 
and modulating the spiral in response to perturbations as 
to maintain a low level of accretio n. The observed l inear, 
Rb ~ 0.6r scaling of bubble sizes (jDiehl et al.l feOOS') may 
provide direct evidence for such mechanical feedback and 
coupling. 

Regardless of the cooling problem and the type of flow, 
a spiral flow strongly affects/mediates the AGN output, 
for example by diverting hot bubbles, distorting jets, and 
deflecting shocks. Indeed, there is evidence that bubbles 
interact with, and appear to be e ntrained by, the spira l 
pattern (see for exa mple Figure [71 iForman etail (|200l : 
iFabian et al.l (|201lD ); such rising bubbles do provide di- 
rect evidence that at least some part of the ICM is an 
outflow. AGN bubbles are therefore natural candidates 
for mediating the mechanical feedback; this can also ex- 
plain the trailing spiral rotation. Interactions with a spi- 
ral flow may explain the rapid g rowth of bubbles and 
their linear size-distance relation ([Diehl et al.ll2008l and 
see £j6j) , t he asymm etric deformation of flow around some 
bubbles (Lim 201 ip. and the large sizes some bubbles at- 
tain (e.g.. iMcNamara et aDl2005f ). 

Bubbles and cold filamentary nebulae, observed near 
the centers of some CCs, gauge the flow. In partic ular, 
coherent axes of bubbles (e.g., Fab ian et al.l 1201 1[) and 
extended linear cold filaments, thought to arise from the 
local thermal instability (e.g., iPizzolato fc Sokerl 120061 : 
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ISharma et al.l 120121) . may seem to limit the presence of 
spiral flows to larger radii. However, (i) most of the vol- 
ume is filled with the slow, < 100 km s _1 gas in which 
differential rotation is minimal, whereas the fast flow is 
confined to a thin layer below the TD; (ii) bubbles appear 
to be entrained by the flow, possibly driving or stabiliz- 
ing it, so extended, albeit slightly bent, axes are to be 
expected; and (Hi) some extended linear structures may 
be parallel to the rotation axis z, for example the long, 
warped filaments sometimes observed in isolation. The 
presence of filamentary nebulae in some CCs may reflect 
enhanced cooling due to the weakening of the spiral out- 
flow in those clusters. 

By focusing on the steady state flow, we have not ad- 
dressed the initial formation of the spiral pattern. How- 
ever, sloshing simulations have demonstrated that such 
patterns easily form in cores in response to strong pertur- 
bations, and standing accretion shock instability (SASI) 
simulations suggest that spiral modes are excited even in 
unperturbed spherical accretion in the presence of cool- 
ing. Note that, if the universal thermal profiles of CCs 
are indeed the signature of a spiral flow, and if these are 
composite inflows (of the first type) of a merger origin, 
then a major merger (of the type invoked to explain the 
minority, ~ 1/3 of non cool-core clusters) would be nec- 
essary in order to generate these spirals. Otherwise, one 
would expect to see clusters with a disrupted core or no 
core at all, with no evidence for a major merger. 

Groups of galaxies show CCs that resemble those in 
clusters, with similar profiles and occasionally CFs, but 
with a higher fraction of gas cooling per unit time (e.g., 
IMcDonald et al.l 1201 ID . Our results suggest that such 
groups too carry spiral flows, but the latter are less effi- 
cient in mixing the gas due to the scaling of AGN feed- 
back with system mass. 

Our analysis bears important implications for inter- 
preting the structure, composition, and gravitational po- 
tential of CCs, showing for example that deviations from 
hydrostatic equilibrium are locally significant. Spiral 
flows may play an important role in mixing the gas on 
large scales, possibly facilit ating the smooth meta llicity 
distribution observed (e.g., ISanderson et al.|[2009f) . and 
the homogen eity of cosmic -ray ions inferred from radio 
observations ([Kcshct 20To|). Our results emphasize that 
spiral flow studies, such as sloshing simulations, should 
incorporate AGN feedback, in particular in the form of 
buoyant bubbles. Such simulations can directly test if 
a spiral flow naturally produces the characteristic core 
structure observed, by studying the evolution and long- 
term stability of flows with different initial conditions. 
Notice that spiral flows cannot be properly simulated 
without accounting for (or at least mimicking the effect 



of) the TD-stabilizing magnetic fields. 

The existence of spiral flows in cores can be directly 
tested through the spectroscopic shift and broadening 
induced by flows along the line of sight. Evidence 
for an intricate, possibly s piral velocity structu re was 
thus uncovered in Perseus ([Sanders et al.l 12004 figure 
8), and evidence for spectral broadening at the level 
of the s ound velocity was f ound in a small subset of 
clusters (| Sanders et al.|[20Tll) . Upper limits on the en- 
ergy associated with the projected velocity dispersion, 
at the level of 5% to 20% of the t hermal energy, were im - 
posed in several othe r clust er s. bvlChurazoy et al.l (2008) ; 
ISanders et all (I20T0L I20II : IBulbul et al. l (|2012D . Thus, 
present spectroscopy is in general unable to probe the 
narrow, fast outflows, whose kinetic energy is a small, 
~ 0.03(r/ro) -8 / 5 fraction of the thermal energy (see 
S|6]). In addition, such observations are only sensitive 
to motions along the line of sight, and typically probe 
scales larger than the base of the spiral, where the flow 
is fastest. 

The spiral bulk flows should be directly observable with 
future telescopes, for example using the next interna - 
tional X-ray satellite, ASTRO-H (jTakahashi et all 20 10:. 
scheduled for launch in 2014. In particular, ASTRO- 
H should be able to spectroscopically resolve the fast 
flows near the centers of nearby, well-resolved clusters, 
in which the spiral is not observed face-on (along z). 

In conclusion, we show that the universal thermal pro- 
files of CCs, their size, multiphase properties, common 
spiral thermal and chemical features, and ubiquitous 
CFs, can all be explained by the presence of bulk spiral 
flows. This suggests that spiral flows may be the funda- 
mental, defining property of CCs. Such flows can glob- 
ally stabilize the core against cooling, while quenching 
mass accretion, by combining a slow, volume-dominant 
inspiral, with a fast, cold, feedback-regulated outflow. 
Present observations are just sensitive to the "tip of the 
iceberg" of this prevalent, large scale phenomenon. 
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